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Obesity is known to be strongly associated with cardio- 
vascular disease and cancer, the leading causes of mortal- 
ity worldwide, and develops owing to interactions between 
genes and the environment. DNA methylation can act as a 
downstream effector of environmental signals, and analy- 
sis of this process therefore holds substantial promise for 
identifying mechanisms through which genetic and envi- 
ronmental factors jointly contribute to disease risk. Global 
DNA methylation of peripheral blood cells has recently 
been proposed as a potential biomarker for disease risk. 
Repetitive element DNA methylation has been shown to be 
associated with prominent obesity-related chronic diseas- 
es, but little is known about its relationship with weight 
status. In this study, we quantified the methylation of Alu 
elements in the peripheral blood DNA of 244 healthy wom- 
en with a range of body mass indexes (BMIs) using py- 
rosequencing technology. Among the study participants, 
certain clinical laboratory parameters, including hemoglo- 
bin, serum glutamic oxaloacetic transaminase, serum glu- 
tamic-pyruvic transaminase, total cholesterol, and triglyc- 
eride levels were found to be strongly associated with BMI. 
Moreover, a U-shaped association between BMI and Alu 
methylation was observed, with the lowest methylation 
levels occurring at BMIs of between 23 and 30 kg/m . How- 
ever, there was no significant association between Alu 
methylation and age, smoking status, or alcohol consump- 
tion. Overall, we identified a differential influence of BMI on 
global DNA methylation in healthy Korean women, indicat- 
ing that BMI-related changes in Alu methylation might play 
a complex role in the etiology and pathogenesis of obesity. 
Further studies are required to elucidate the mechanisms 
underlying this relationship. 
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INTRODUCTION 

DNA methylation is a core epigenetic process that influences a 
wide variety of biological mechanisms, including gene expres- 
sion, chromosomal stability, imprinting, and cellular differentia- 
tion (Bernstein et al., 2007). Methylation can mediate environ- 
mental influences on gene expression and can modulate dis- 
ease risk associated with genetic variation (Foley et al., 2009). 
Abnormal DNA methylation patterns (including genome-wide 
hypomethylation and gene-specific hypo- and hypermethyla- 
tion) have been shown to be associated with a range of health 
outcomes (Ozanne and Constancia, 2007). In cancer and cer- 
tain other diseases, most of these changes have been ob- 
served at the tissue level (Esteller et al., 2001 ; Rodenhister and 
Mann, 2006). Data relating to whether DNA methylation chang- 
es in peripheral blood cells can serve as useful, informative 
biomarkers for different health outcomes is much more limited, 
but new information is rapidly emerging (Heyn and Esteller, 
2012). A number of methods are available for the analysis of 
global DNA methylation levels [total content of 5-methyl- 
cytosine (5-mC)] (Terry et al., 2011). There are over 500,000 
long-interspersed nuclear elements (LINEs) and 1,500,000 
short-interspersed nuclear elements (SINEs) across the human 
genome. Alu elements are the most common SINE; comprising 
11% the human genome and containing a third of all methyla- 
tion sites (Deininger and Batzer, 1999). The substantial repre- 
sentation of these repetitive elements throughout the genome 
makes their effectiveness as proxies for global DNA methyla- 
tion highly evident. Indeed, methylation levels in Alu and LINE-1 
repeats have been shown in previous studies to be associated 
with total genomic methylation content (Weisenberger et al., 
2005; Yang et al., 2004). Estimates of these are increasingly 
being used in epidemiological investigations owing to their rela- 
tively low-cost, high-throughput, and the fact that they can offer 
quantitative results (Weisenberger et al., 2005). Interestingly, 
global methylation measures in blood DNA have been shown to 
vary in relation to demographic and lifestyle characteristics, 
including age, sex, smoking status, alcohol consumption, phys- 
ical activity, and diet (Heyn and Esteller, 2012; Zhu et al., 2012). 
However, there are still few studies of the relationship between 
DNA methylation and weight status. 

Obesity is simply defined as a condition of abnormal or ex- 
cessive fat accumulation in adipose tissue, to the extent that 
health may be impaired. Obesity is the result of the interplay 
between external (environmental) and internal (genetic) factors. 
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It has been cited as the epidemic of our time, with obesity rates 
sharply and steadily rising in many parts of the world, leading to 
increased morbidity and mortality due to type 2 diabetes and 
cardiovascular disease (Danaei et al., 2009). Obesity-related 
mortality by far exceeds mortality caused by other common 
diseases (Anderson and Caswell, 2009). Moreover, about 20% 
of all cancers are caused by overweight, and obesity has been 
shown to be associated with increased risk and worse out- 
comes following diagnosis (Wolin et al., 2010). Furthermore, 
the incidence of some types of cancer have been shown to be 
strongly associated with body mass index (BMI), which has 
historically been used to define obesity (Renehan et al., 2008). 
Our understanding of how and why obesity develops is incom- 
plete. To date, several studies concerning DNA methylation and 
obesity have focused primarily on gene-specific methylation 
and had relatively small sample sizes (Carless et al., 2013; 
Feinberg et al., 2010; Fujiki et al., 2009; Milagro et al., 2012; 
Wang et al., 2010). Currently, studies of global methylation 
levels in blood DNA and BMI have yielded inconsistent results 
(Kim et al., 2009; Perng et al., 2013; Piyathilake et al., 2011; 
Wang et al., 2010; Zhang et al., 2011a;2011b; Zhu et al., 2012). 
BMI is now considered to be an important determinant of meth- 
ylation biomarkers in the blood of women of reproductive age 
(van Driel et al., 2009). Interestingly, recent evidence has 
shown that alterations in global DNA methylation may be an 
important contributor to incidence risk of cancer and cardiovas- 
cular disease (Baccarelli et al., 2010; Terry et al., 2011). In this 
study, in order to investigate the specific effects of BMI on glob- 
al DNA methylation and develop blood-based biomarker as 
initial screening test for obesity-related disease, we quantified 
the methylation of Alu elements in the peripheral blood DNA of 
244 women with a range of BMIs using pyrosequencing tech- 
nology. 

MATERIALS AND METHODS 

Study population 

Study subjects were 244 apparently healthy women aged be- 
tween 20 and 51 years (mean 32 ± 7.9 years). This study was 
approved by the Institutional Review Board of Kyungpook Na- 
tional University Hospital. Additionally, informed written consent 
was obtained from all subjects before they participated in the 
study. Demographic information and lifestyle factors were de- 
termined for all participants by trained interviewers using a 
standardized questionnaire via face-to-face interviews. Height 
and bodyweight were measured using standard methods with 
participants wearing light clothes. BMI is easy to obtain [calcu- 
lated by weight divided by height squared (kg/m 2 )] and is con- 
sidered to be an acceptable proxy for under and overweight. 
BMI has also been shown to be directly related to health risks 
and mortality in many populations. According to the current 
international standard (WHO 1998), a BMI of over 25 kg/m 2 is 
considered to indicate overweight and a BMI of over and 30 
kg/m 2 is considered to indicate obesity. However, the BMI cutoff 
values for overweight and obesity vary in Asian populations 
from 22 to 25 kg/m 2 and the increased risk of co-morbidities 
conferred by obesity have been shown to occur at a lower BMI 
in several Asian populations (Low et al., 2009). We therefore 
divided the participants into 3 categories based on BMI normal 
weight (BMI < 23 kg/m 2 ), overweight (23 kg/m 2 < BMI < 30 
kg/m 2 ), and obese (BMI > 30 kg/m 2 ). Blood samples were ob- 
tained via venipuncture after overnight fasting, and serum sam- 
ples were separated by centrifugation and transferred to uncon- 
taminated bottles with Teflon-coated caps. All samples were 



kept frozen at -70°C until analyses were conducted. Clinical 
laboratory values were determined by standard biochemical 
automatic or semi-automatic methods. 

Blood DNA extraction and bisulfite treatment 
Genomic DNA was extracted from whole-blood samples using 
a QIAamp DNA Blood Kit (Qiagen, USA). One microgram DNA 
was bisulfite-modified using an EZ DNA Methylation-Gold Kit 
(Zymo Research, USA) according to the manufacturer's in- 
structions. Final elution was performed with 30 ul M-Elution 
Buffer (Zymo Research) and was stored at -70°C until analyzed. 
Built-in analysis of non-CpG cytosine residues provided an 
internal control for the completeness of bisulfate treatment. 
Following treatment with bisulfite, the conversion of this C into T 
is expected to be 100%. It is possible to insert a C/T single- 
nucleotide polymorphism into the sequence to be analyzed, 
and this will result in 100% Ts if conversion is efficient. 

Assay of Alu methylation 

Alu methylation analysis was quantitatively performed on the 
bisulfite-treated DNA using pyrosequencing with the polymer- 
ase chain reaction (PCR) primers and conditions previously 
described (Bollati et al., 2007). In brief, the bisulfite-treated 
samples (50 ng) were amplified with a biotin-labeled primer via 
PCR, which enables the conversion of the PCR product to a 
single-stranded DNA template suitable for pyrosequencing. 
Confirmation of the quality of the PCR products and their free- 
dom from contamination was established on 2% agarose gels 
with ethidium bromide staining. After purification of PCR prod- 
ucts using Sepharose beads on PyroMark Vacuum Prep Work- 
station (Qiagen), pyrosequencing was carried out using the 
PyroMark Q96MD System (Qiagen) according to manufactur- 
er's instructions, including a single-strand binding protein. The 
percentage of methylatation was expressed for each DNA locus 
as % 5-mC divided by the sum of methylated and unmethylated 
cytosines. We tested each marker 3 times and used the aver- 
age percentage in the statistical analysis. Three controls were 
included in every pyrosequencing run. One well was filled with 
water to ensure no contamination, and 2 wells were filled with 
CpGenome Universal methylated and unmethylated DNA 
(Chemicon, USA) to evaluate the repeatability of the assay. 

Statistical analysis 

Statistical analysis and plotting was performed using R version 
2.15.3 (http://www.r-project.org). Continuous variables were 
compared using one-way analysis of variance or analysis of 
covariance and pairwise post hoc comparison identified statisti- 
cal differences following Bonferroni adjustment. Categorical 
variables were compared using a chi-squared test. 

RESULTS AND DISCUSSION 

Growing evidence indicates that there are sex differences in 
global DNA methylation in the blood (Terry et al., 2011; Zhu et 
al., 2012). We have chosen healthy women as study subjects. 
First, overweight and obesity pose more serious health chal- 
lenge to women than men. Obesity is a common disorder af- 
fecting approximately 1 in 3 women and the prevalence of ex- 
treme obesity is about 50% higher among women than men 
(Ryan and Braverman-Panza, 2014). Moreover, prevalence of 
metabolic syndrome is associated with obesity and its escala- 
tion has been steeper in women, particularly in young women 
(Regitz-Zagrosek et al., 2007). Second, recent result has shown 
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Table 1. Characteristics of study participants 





Normal weight 


Overweight 


Obese 


F/)C 


P-value 


Post hoc compar- 




(n = 92) 


(n = 79) 


(n = 73) 






ison using the 

LJUilGdM lllcHMOU 


Age, years (SD) 


31.52 (8.84) 


31.70 (5.52) 


32.93 (8.98) 


0.72 


0.487 




Smoking status, n (%) 
Non-smoker 
Former smoker 
Current smoker 


R7 170 a\ 

VI \l £..0) 

6 (6.5) 
19(20.7) 


60 (75 
6 (7.6) 
13(16.5) 


52 (71 2) 
2 (2.7) 
19(26.1) 


3.51 


0.477 




Alcohol consumption, n (%) 
Non-drinker 
Ever drinker 


22 (23.9) 
70 (76.1) 


23(29.1) 
56 (70.9) 


30 (41.1) 
43 (58.9) 


7.58 


0.270 




Hemoglobin, % (SD) 


11.97(1.05) 


12.23 (0.99) 


12.57(1.22) 


6.18 


0.002 


3 > 2, 1 


Creatinine, mg/dl (SD) 


0.87 (0.12) 


0.86 (0.09) 


0.87 (0.10) 


0.53 


0.587 




Glucose, mg/dl (SD) 


93.72 (6.17) 


95.57(23.19) 


93.18(14.18) 


0.49 


0.611 




GOT, units/I (SD) 


17.50(4.48) 


19.67 (7.17) 


24.27(16.43) 


9.08 


< 0.001 


3 > 2, 1 


GPT, units/I (SD) 


13.36 (8.08) 


18.48(12.58) 


29.92 (27.62) 


18.80 


<0.001 


3 > 2, 1 


Total bilirubin, mg/dl (SD) 


0.62 (0.26) 


0.59 (0.30) 


0.58 (0.27) 


0.71 


0.492 




Total cholesterol, mg/dl (SD) 


179.59 (31.66) 


190.94(47.89) 


199.08 (30.35) 


5.68 


0.004 


3 > 2 > 1 


Triglyceride, mg/dl (SD) 


83.93 (46.07) 


116.57 (86.88) 


168.48(128.18) 


17.92 


< 0.001 


3 > 2> 1 



SD, standard deviation; GOT, glutamic oxaloacetic transaminase; GPT, glutamic pyruvic transaminase 




Fig. 1. Scatter plot and simple regression of 
population characteristics and body mass 
indexes (BMIs). (A) hemoglobin, (B) serum 
glutamic oxaloacetic transaminase (GOT), 
(C) serum glutamic pyruvic transaminase 
(GPT), (D) total cholesterol, (E) triglyceride 
levels. 



that BMI is a strong determinant of methylation biomarkers in 
women of child-bearing age (van Driel et al., 2009), indicating a 
greater likelihood of BMI effects on global DNA methylation in 
women. Finally, previous studies investigating the effect of BMI 
on global DNA methylation include mainly male populations 
(Zhu et al., 2012). In addition, we have found that Alu assays 
may be a more sensitive method for detecting global DNA hy- 
pomethylation in peripheral blood leukocytes in the general 
population than LINE-1 assays (Kim et al., 2009). Most studies 
have focused primarily on LINE-1 methylation as a surrogate 
for global DNA methylation (Perng et al., 201 3; Piyathilake et al., 
2011; Zhang et al., 2011a; 2011b). In this study, we determined 



Alu methylation levels in a population of Korean women. The 
characteristics of the study population are shown in Table 1 . 
Obese participants had higher baseline hemoglobin, serum 
glutamic-oxaloacetic transaminase (GOT), serum glutamic- 
pyruvic transaminase (GPT), total cholesterol (TC), and triglyc- 
eride (TG) levels than normal weight and overweight partici- 
pants. Scatter plot and linear regression analyses showed that 
these biochemical parameters were positively correlated with 
BMI (Fig. 1), which is consistent with recent observations that 
GOT and GPT levels are strongly associated with BMI, TC, and 
TG in Asian populations (Lin et al., 2010; Sull et al., 2009). 
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Table 2. Correlation between Alu methylation and BMI 



Gene 




Crude means 






Adjusted means* 




Mean 


F 


P-value 


Mean 


F 


P-value 


Alu methylation 














Normal weight 


26.27 + 0.82 


67.46 


< 0.001 


26.28 + 0.09 


67.72 


< 0.001 


Overweight 


24.96 + 0.45 






24.95 ±0.12 






Obese 


25.99 ± 0.94 






25.96 ±0.10 







*Adjusted for clinical laboratory values including hemoglobin, serum glutamic oxaloacetic transaminase, serum glutamic pyruvic transaminase, total 
cholesterol, and triglyceride levels. The Bonferroni method was used for post-comparison. 



B 




• BMI < 23 0 
o 23 0*BMI<30 0 
■ 30 0 s BMI 



Fig. 2. Graphic representation of body mass 
index (BMI)-related global DNA methylation. 

(A) Scatter plot of BMI and Alu methylation, 

(B) Percentage of Alu methylation of blood 
cells in normal weight, overweight/pre-obese, 
and obese subjects. 



Normal Weight Oveiweighl or Preobese 



Taken together, our data support the hypothesis that body- 
weight may be the major factor in determining serum liver en- 
zyme levels. Moreover, a significantly lower level of Alu methyl- 
ation was found in overweight participants compared with nor- 
mal weight or obese participants (p < 0.0001 ), but no significant 
difference was observed between the latter 2 groups, indicating 
a U-shaped association between Alu methylation and BMI (Ta- 
ble 2 and Fig. 2). After adjusting for other clinical covariates, this 
association was still found to be significant (Table 2). However, 
age, smoking status, and alcohol consumption were not signifi- 
cantly associated with Alu methylation levels in blood cells (data 
not shown). Similarly, previous studies have reported no asso- 
ciation between smoking, alcohol consumption, and leukocyte 
DNA methylation (Hsiung et al., 2007; Kim et al., 2009; Terry et 
al., 2011; Zhuetal., 2012). 

To date, studies of global methylation levels of blood DNA 
and BMI have yielded inconsistent results (Kim et al., 2009; 
Pemg et al., 2013; Piyathilake et al., 2011; Wang et al., 2010; 
Zhang et al., 2011a; 2011b; Zhu et al., 2012). Most previous 
studies reported no significant association between BMI and 
global DNA methylation (Kim et al., 2009; Wang et al., 2010; 
Zhang etal., 2011a; 2011b; Zhu etal., 2012). In contrast, higher 
BMI has been reported to be associated with lower LINE-1 
methylation (Pemg et al., 2013; Piyathilake et al., 2011). Nota- 
bly, we observed that participants with a BMI of 25 to 30 kg/m 2 
had the lowest methylation levels compared with the other 
groups, and that Alu methylation was elevated among those 
with a BMI higher or lower than that range. To the best of our 
knowledge, this is the first study to report a novel U-shaped 
association between BMI and Alu methylation. Interestingly, the 
association between BMI and all-cause mortality has previously 
been shown to be U-shaped, with the concave region corre- 
sponding to BMIs of between 22 and 26 kg/m 2 (Berrington de 



Gonzalez et al., 2010; Whitlock et al., 2009; Zheng et al., 2011). 
Recently, we have observed that persistent organic pollutants, 
polychlorinated biphenyls are related to MGMT hypermethyla- 
tion with an inverted U-shaped curve (P for quadratic term < 
0.01). The prevalence of hypermethylation was highest in the 
subjects in the 2 nd quintile (28.4%); however the subjects in the 
5 th quintile had a low prevalence of hypermethylation (8.2%) 
(unpublished data). Taken together, these data suggest the 
biphasic (hermetic) dose response model in which adaptive 
responses to low doses otherwise harmful conditions improve 
the functional ability of the cells and organism, resulting in a U 
or inverted U-shaped dose response (Calabrese, 2010). Ac- 
cordingly, several lines of evidence are presented that DNA 
methylation may be an attractive candidate molecular mecha- 
nism for hormesis-like response (Vaiserman, 2011). 

Although the exact mechanism underlying the elevated Alu 
methylation observed in the blood cells of obese people is not 
known, several possible explanations have been suggested. 
First, it is reasonable to speculate that sleeping disturbance 
might contribute to increases in Alu methylation in obese peo- 
ple. Obstructive sleep apnea (OSA) is highly prevalent among 
obese patients (Young et al., 1993) and methylation of the en- 
dothelial nitric oxide synthase and forkhead box P3 genes in 
blood DNA has been shown to be linearly associated with the 
severity of OSA (Kheirandish-Gozal et al., 2013; Kim et al., 
2012). This suggests that global methylation can be influenced 
by OSA and other sleep disorders. In support of this theory, we 
found that scored sleep management was higher in normal and 
overweight subjects compared with obese ones (data not 
known), indicating that obese individuals may have poorer 
sleep quality than non-obese people. Second, although BMI 
does not differentiate between muscle mass and adipose tissue, 
people with high BMI can be inferred to have high levels of 
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body fat, causing higher levels of inflammatory markers such as 
interleukin (IL)-6 (Ramkumar et al., 2004). Moreover, persistent 
inflammation may cause global DNA hypermethylation in pe- 
ripheral blood cells via IL-6 signaling (Stenvinkel et al., 2007). 
Thus, it is likely that elevated Alu methylation in people with 
high BMI is a consequence of a low but continuous inflammato- 
ry drive of obesity. Third, Kamei et al. (2010) demonstrated that 
the expression of DNA methyltransferase 3a is markedly up- 
regulated in the adipose tissue of obese mice. Similarly, a sig- 
nificantly higher frequency of RASSF1A, p16 and hMLH1 
methylation is observed in nonmalignant tissues of obese pa- 
tients (Peters et al., 2007; Ye et al., 2006). In light of evidence 
that tissue-specific methylation profiles are reflected in circulat- 
ing cell epigenomics (Yang et al., 2010), it is tempting to specu- 
late that obesity can increase global methylation levels in pe- 
ripheral blood DNA. Fourth, there is compelling evidence that 
heterogeneity in the number of specific leukocytes has the 
potential to confound methylation measurements using whole 
blood DNA and that increased total leukocytes are associated 
with obesity (Talens et al., 2010; Veronelli et al., 2004). Thus, a 
plausible hypothesis is that a progressively uprising trend at 
high BMI might be due to differential lymphocyte counts associ- 
ated with adiposity development. In this respect, it is noteworthy 
that DNA methylation is closely related to lymphocyte count 
(Zhu et al., 2012). Finally, these divergent results may be due to 
the specific racial/ethnic characteristics of the study population 
or the small sample size. Studies with a larger sample size are 
therefore required to confirm these findings. Further studies to 
confirm the U-shaped association between Alu methylation and 
BMI by measuring LINE-1 methylation and to understand the 
biological mechanisms underlying observed relationship are 
also required. 

Collectively, the present study showed a novel U-shaped as- 
sociation between Alu methylation and BMI, suggesting that 
BMI-global DNA methylation relationship might play a convolut- 
ed role in etiology and pathogenesis of obesity as a hermetic 
dose response. In addition, the correlation could be a potential 
biomarker for the initial screening of other obesity-related dis- 
ease. Furthermore, our results showed magnitude of obesity in 
general healthy population as Korean reproductive aged wom- 
en. 
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